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Abstract 


The helium abundance in the Jovian atmosphere is derived from Voyager 1 
data by two methods. The first method uses only infrared spectra from 
selected locations on the planet while the second method uses a thermal 
profile independently derived from radio occultation measurements and infrared 
spectra recorded near the jccultation point. A hydrogen mole fraction of 
0.897 t 0.030 is obtained ^rom the first method while the second method gives 
0.880 ± 0 . 036 , corresponding to helium mass fractions of 0.19 t 0.05 and 0.21 
t 0.06 respectively. The estimated errors for the* first method are primarily 
due to systematic uncertainties in the and He absorption coefficients, 
while those for the second method result mainly from errors in the radio 
occultation profile and are lees well known. Random errors in the measured 
infrared spectra are found to be negligible in both cases. The results are 
consistent with a uniform mix of hydrogen and helium within Jupiter's 
interior, but a modest amount of helium depletion (aY < 0.05) cannot be 
excluded. 
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Introduction 


An accurate determination of the helium abundance in the atmosphere of 
Jupiter is of great importance to models^ of the planetary interior and of the 
solar system evolution. Jupiter's large planetary mass, combined with a low 
exospheric temperature, has prevented significant escape of hydrogen and 
helium since planetary formation, and the bulk composition of Jupiter should 
be similar to that of the primitive solar nebula. Therefore, a oompar son of 
the atmospheric abundance of helium on Jupiter with estimates of solar and 
cosmic abundances may Indicate » tether helium differentiation has occurred 
within the planet. 

Pressure induced absorption from Hg-Hg and Hj-He collisions accounts for 
a large fraction of the far infrared opacity of Jupiter (Trafton, 1967). 

Taking advantage of the different spectral dependences of the H 2 -H 2 and H 2 ~He 
absorption coefficients, Gautier and Grossman (1972) proposed a method of 
Inferring the helium abundance from spectral measuremv “s in the 300-700 cm“^ 
region. This method has been used to estimate the hydrogen mole fraction for 
Jupiter from airborne-spectroscopic (Houck, et al., 1975; Gautier, et al., 
1977) and Pioneer center-to-limb measurements (Orton and Ingersoll, 1976). 
However, the accuracy of the results was limited in the airborne measurements 
by the absolute calibration and lack of spatial resolution and in the Pioneer 
determination by the use of broad band filters and the Implicit assumption of 
horizontal homigeneity in the temperature field. The Voyager IRIS data provide 
significant advantages over previous measurements. A large number of 
calibrated spectra were obtained from various locations on Jupiter with a 
spectral resolution of 4.3 cm”^. This permits the selection of cata from 
locations where the effect of aerosols is expected to be small. 

In a preliminary analysis of Voyager data, Hanel et al. (1979) estimated 
the helium mole fraction to lie between 0.08 and 0.14. In the present paper 
the estimate is refined through more extensive analyses following two 
different approaches. First, the inversion method of Gautier and Grossman 
( 1972 ) is applied to the IRIS infrared spectra from clear areas; next, radio 
occultation profiles are combined with IRIS data acquired near the ocoultation 
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points. Tht sources of error associated with each approach are examined, and 
the implications of the results are 


Collision-Induced Absorption of 

Both methods used in this investigation depend o'i calculations of 
synthetic Jovian spectra in the 280-600 cm**^ region. Tlierefore, before 
proceeding to a discussion of the methods themselves, the formulation of 
collision induced H 2 absorption coefficients used in these calculations will 
be reviewed. The formulation is based on the recent laboratory and 
theoretical work of Birnbaum and Cohen (1?76), Birnbaum (1978), and Cohen and 
Birnbaum (1980). 

The H 2 absorption coefficients enter the radiative transfer equation 
through the expression for monochromatic optical depth. 


s ;[q^ A ♦ q(1-q) B] (P/P(jxp)^ 


where q is the nole fraction of H 2 , P is the total atmospheric pressure, T is 
the atmospheric temperature, and z is the geometric altitude. The hydrogen 
mole fraction is defined by 


q 




n 


He 


where n is the number density of the corresponding gases. The collision- 
induced absorption coefficients for H 2 -H 2 and H 2 “He collisions are represented 
by A and B, respectively. 


Molecular hydrogen has a center of symmetry and, consequently, no 
permanent dipole moment. However, a weak collision-induced dipole spectrum 
exists which creates significant absorption for the very long path lengths 
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•noount«r«d in the Jovltn atmosphere; about 40 km atm are above the one bar 
level. The theory of oollislon-induced absorption in homonuclear molecules 
indicates the existence of an induced dipole which results from two distinct 
physical mechanisms. The first mechanism, qusdrupole interaction, is 
relatively long range and anlsotronic in nature. The second mechanism, 
electron overlap interaction, is relatively short range and has both isotropic 
and anisotropic terms. The major contribution to the absorption for H 2 -H 2 
collisions is from the quadrupolar interaction, with a smaller contribution 
from anisotropic overlap. The isotropic overlap is forbidden by symmetry for 
H„-H- collisions, but for H--He collisions it is the dominant factor. The 
quadrupolar and anisotropic overlap terms a 'e approximately one tenth of the 
terms tor H 2 -H 2 collisions because of the smaller polarizability of He. 

The analytical collision-induced llneshape of Birnbaum and Cohen (1976) 
was used for A and B, along with the appropriate shape factors determined from 
laboratory oata (Cohen and Birnbaum, 1980). The expression for A is 

A oc V l‘o [# r(v, r“) 

+ r”) 


with V being the transition wavenumber, the summed intensity of the 
translational-rotational transitions due to contributions from the quadrupole 
and anisotropic overlap induced dipole, p a Boltzmann I'actor, 4 > a sum over 
states, r a lineshape function, and and T 2 lineshape parameters controlling 
the width and shape of r. The 11 superscript is for H 2 -H 2 collisions. The 
first term in the brackets represents absorption in the translational modes, 
caused by a change in relative kinetic energy of two colliding H 2 molecules, 
while the remaining terms represent absorption in the rotational-translational 
modes by the resonance lines. All H 2 energy levels are assumed to be in local 
thermodynamic equilibrium. 
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The A absorpcion coefficient is illustrated in Fig. la and b for 
temperatures of 180 and 100 K, respectively. This range approximates the 
temperatures of atomospheric layers in which the Hp lines are formed. The 
effect of temperature on the relative contribution from the transnational mode 
(0-20C cm"^), the S(0) resonance line (354 cm”^) and the S(1) resonance line 
(587.7 CiTi"^') may be seen. At these low temperatures only the 2 lowest 
rotational states, J > O (para) and J » 1 (ortho), are significantly populated. 

The B absorption coerficient for H^-He collisions is represented by 


where the superscript 12 refers to H^-He collisions and the terms are as 
defined before. The second term in the rectangular bracket represents the 
absorption due to the isotropic overlap with being the overlap intensity. 

The B absorption coefficient is also illustrated in Figs, la and b. The 
weaker absorption of the S lines, relative to A, is due to the weaker 
quadrupole induced dipole for He whi^e the stronger translational band arises 
from the isotropic overlap term in B. This term leads to the enhancement of 
the Hg-Hp absorption coefficient which allows the helium abundance to be 
derived from the infrared spectrum. 

The influence of He on the optical depth of H 2 in the Jovian atmosphere 
may be demonstrated by the ratio 

R = T(q)/T(q') = 0.9 

where the value q* s 0.9 is taken as representative of the true Jovian 
hydrogen mole fraction. This ratio is shown in Fig. 2 for a mean atmospheric 
temperature of 140 K, and for q = 1.0 and 0.8. For wavenumbers greater than •/' 
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300 om'^ R shows a strong dependence on q, but only a weak dependence on 
wavenumber. This behavior refelcts the dominance of the ^ 2~^2 
expression for the optical depth. If the atmospheric temperature is known 
independently as a function of the total atmospheric pressure, the measured 
spe'^tral radiance can be regarded as a function of q alone, and this strong 
sensitivity can then be exploited. This forms the basi.s of the combined 
IRIS-radio ocoultatlon method. The stronger wavenumber dependence of R in the 
200 cm"' region refelcts the increased relative importance of the Hj-He term 
in the optical depth. This frequency dependent sensitivity to q forms the 
basis of the inversion method in which both the atmospheric temperature 
profile and q are retrieved from the spectral radiances dlone. Unfortunately, 
the region of strongest wavenumber dependence is not accessible in the Jovian 
infrared spectrum because of Interference from the rotational NH, lines. In 
the accessible region 280-600 cm , the wavenumber dependence is relatively 
small; however, with the large number and good quality of spectra obtained by 
the Voyager IRIS it is still possible to extract significant information as 
discussed in the following section. 


Inversion Metnod 


The principle of the Inversion method as applied to Jupiter can be 
understood from the following considerations. Within the S(0) and S(1) lines 
of the spectrum generally two or more wavenumbers exist for which the 
brightness temperatures are the same; see for example Fig. 8. Therefore, 
measurements at these wavenumbers must be redundant in the sense that unit 
optical depth occurs at approximately the same atmospheric pressure level, and 
the same portion of the temperature profile is sampled in each case. However, 
a small but detectable difference in the dependence of optical depth on q 
exsists for each wavenumber as discussed In the preceeding section. For this 
reason, the dependence on atmospheric temperature can be eliminated; only one 
value of q satisfies a given set of measurements at the redundant wavenumbers. 
In practice, temperature profiles are retrieved for various values of q from 
sets of radiances which include redundant points, for which the corresponding 
weighting functions are nearly coincident. The solution for q is then taken 
as that value for which the rms residual between measured and calculated 
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radiances is a minimum, The wavenumbers chosen for this analysis are 267# 

310, 3«0, 365. U75, 530, and 602 cm‘^. 

Unrecognized sources of opacity such as clouds may cause errors. Since 
regions of strong 5 urn emission are known to be relatively cloud free, only 
spectra with brightness temperatures greater than 250 K at 5 urn were Included 
In the analysis. On the basis of this criterion, 62 spectra were selected, 
the majority of which were from the North Equatorial Belt. 

Each of the chosen spectra was Inverted for various values of q, and the 
root mean square of the residuals was calculated as a function of q. A 
typical result Is shown In Fig. 3 where the minimum occurs at q > 0.89, and a 
histogram of all minima of q for the 62 samples Is given In Fig. U. From this 
distribution we obtain a mean value of q i 0.897 with a standard deviation a = 
0.015 and a standard error of the mean Aqj^ = 0.002. The small value of Aqj^ 
Indicates that the random errors due to the Instrument noi;,8 have a negligible 
effect. However, systematic errors can be present due to calibration 
uncertainties and errors In the absorption coefficients. 

The accuracy of the IRIS calibration depends on the stability of the 
instrument temperature (Hanel et al., 1981). During the time of the 
measurtnents the interferometer temperature was stable and the uncertainty in 
the instrument temperature was less than 0.1K. Such an effect is negligible 
for the retrieval of q. 

According to Blrnbaum (1978) errors in the absorption coefficients 
probably do not exceed ±5%. The most Important effects for this work are 
attributable to systematic ei ••ors In the H 2 -H 2 and H 2 “He components of the 
coefficients. Several possible forms of error have been examined. A 
frequency independent variation of both A and D of ± 5X, is found to result in 
an error Aq of t 0.005. The introduction of a frequency independent relative 
error in A or B alone produces a much stronger effect. Fig. 5 shows the error 
Aq as a function of the relative error aB/B; a similar diagram will result 
from aA/A. Assuming upper limits for aA/A and aB/B of i 0.05 we obtain an 
uncertainty for Aq of t 0.03. 
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Our final •stimata of th« hydrogtn molar fraction from tha Invarslon 
mathod Is q • 0.697 i 0.030. Aa dlacussad, tha unoartalnty Is primarily due 
to oonsa*^vatlvaly astlmatad systamatlo arrors In tha absorption ooafflclants, 
a factor oftan naglaotad by othar authors. 


IR IS-Radlo Occultatlon Mathod 


Tha Voyagar Radio Sclanca team has derived temperature profiles from 
measurements acquired during occultatlon of thn spacecraft by the Jovian 
atmosphere (Eshleman et al., 1979). Tha upper atmosphere temperatures 
obtained with this technique depend strongly on the boundary condition assumed 
at the top of the atmosphere, but below the tropopause tne solutions become 
increasingly insensitive to the assumed boundary condition. Thirteen IRIS 
spectra were obtained in the vicinity of the Voyager 1 ingress occultatlon 
point (9° to 12°S, 59° to 6?°W) with a field of view of 16° in latitude and 
longitude projected on the planet. This area is considerably larger than the 
narrow channel the radio signals propagate through. The thermal profile 
retrieved from the infrared data by means of inversion techniques described by 
Conrath and Gautier (I960) is shown in Fig. 6. Also shown are temperature 
profiles derived from the radio occultatlon measurements (G. Llndal, private 
communication) for two different boundary conditions at the top of the 
atmosphere. The agreement between the occultatlon profiles and the IRIS 
profiles is good in the troposphere and acceptable in the stratosphere when 
the limited vertical resolution of the Infrared measurements in this portion 
of the atmosphere is taken into consideration (Gautier and Revah, 1975; Manei 
et al., 1977; Orton, 1977). The tropospheric agreement also suggests that the 
presence of clouds and haze dio not have a strong effect on the profile 
derived from the Infrared spectra. 

The radio occultatlon profiles are derived under the assumption of a 
specific hydrogen mole fraction. Changing this value rescales the temperature 
and pressure at each level according to the relations 


and 


P ^ 

^ const 

Warm ocoultatlon profiles derived for q ■ 0.94 and 0.06 are shown in Fig. 7. 
This behavior provides a sensitive method for determination of the hydrogen to 
helium ratio. Radio ocoultatlon profiles are generated for several values of 
q and each profile is then used to cslculate a corresponding spectrum which is 
compared to the measured spectrum. A vslue of q corresponding to the best fit 
is then selected. The calculated spectra are hif nly sensitive to the hydrogen 
mole fraction because a change in q results in changes in both the temperature 
profile and atmospheric transmission which set in the some sense. Only a weak 
dependence exists on the upper atmospheric boundary since the emission in the 
285-600 cm”^ spectral region originates predominantly in the troposphere. 

Determinations of q h9^ e b?en made using both the ingress and egress 
profiles as given by Linca- , et al. (1981). To demonstrate the sensitivity of 
the technique, results of calculations .'or the Ingress profile are shown in 
Fig. 6 where the best fit is obtained fo«* q > 0.875. Similar calculations for 
the egress profile were compared with an average of 14 IRIS spetra acquired 
near the egress point (0.7°N, 313°W), and the best fit was obtained for q s 
0.880. 


The values of q obtained by this method are subject to severe' sources of 
error, including noise in the Infrared spectra, calibration errors, inaccurate 
knowledge of absorption coefficients, and uncertainties In the occultatlon 
temperature profiles as well as possible interference from haze and clouds. 
Each of these sources will now be considered. 

The noise equivalent spectral radiance (NESR) of Voyager 1 IRIS in the 

-1 -9 -2 -1 -1 

280-600 cm spectral region does not exceed 5.i0 watts cm sr /cm 

(Hanel et al., 1980) which results in a signal to noise ratio per spectrum of 

about 300 at 285 cm”^ and 30 at 600 cm”'. These values correspond to errors 

in brightness temperatures of about 0.15K at 285 cm” and 0.5K at 600 cm” . 

An uncertainty of IK in the calculated or measured brightness temperature 

results In an error of about 0.01 in q; therefore, if we make the conservative 
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assumption that ths srrors assoclatad with an Individual spaolrum also apply 
to tha 13 avaragad spaotra, an uppar limit on tha unoartainty in q dua to tha 
random noisa is aq ■ i O.OOU. 

Durin »5 tha maasuramant of tha oocultation points, tha inter faromater 
tamparatura was axtramaly stabla, and accordingly tha absolute calibration 
errors are astimatad to ba lass than tha random errors and have a negligible 
affect on tha determination of q. 

To_ estimate the effect of absorption coefficient errors, outgoing spectra 
have bean calculated with tha total absorotion coefficient modified by t5A>. 

The results in terms of brightness temperature are shown in Fig. 9. The 
amplitude of the effect depends on tha lapse rata in the atmospheric region 
where the emission originates. The effect is weak jt 600 cm"\ which 
corresponds to the tropopause level and is a maximum at 270 cm”^ which 
corresponds to the beginning of the convective zone. These brightness 
temperature changes imply that absorption coefficient errors result in an 
uncertainty in q of about 4q * ± 0.005. 

Tne effect of hazes and clouds may be noti';atle at spectral intervals in 

which the emission emerges from deep atmospheric layers, such as the 450-500 
-1 

cm range and below 300 cm . This may account for the residual differences 
over •’hese intervals between calculated and measured spectra in Fig. 8. 

However neither haze nor clouds are expected to significantly affect the broad 
centers of the 3(0) and 3(1) lines near 350 and 600 cm“\ Therefore the match 
between measured and calculated spectra was optimized for the spectral regions 
tteai’ the line centers. Thus, the errors due to unrecognized haze and cloud 
effects are expected to be small. However, no quantitative estimate has been 
made. 


Tne principal uncertainty in the determination of q comes from the 
uncertainty in the occultation profiles. Lindal, et al. (1981) indicate that 
the tropospheric portion of the ingress profile is leso reliable than the 
egress profile because of the presence of strong turbulence below the 300 mbar 
level at the ingress occultation point. For this reason only the egress 



profil* will b« oonsid«rtd furthtr In th4 present analysis. Estimation of ti'e 
uncertainties In the egress profile la difficult since a detailed error 
analysis has not been published. Llndal et al. quote an error of ± 7K at 1000 
mbar and Indicate the main error source la uncertainty In q. Fo purposes of 
this analysis we shall assume that In the absence of an uncertainty In q, the 
error In temperature would be i 3.5K which la half of the quoted total error. 

This value Is also consistent with an error estimate of 5 to 10K at the 3o 
level given previously by Eshelman, et al. (1979). An error of t 3.5K In the 
oooultatlon profile results In an uncertainty of t 0.035 In the Inferred value 
of q. 


If the various uncertainties discussed above are treated as statistically 
Independent sources of error, then the final result obtained with this method 
Is q » 0.880 i 0.036. 

Discussion and Conclusions 


The hydrogen mole fraction for the atmosphere of Jupiter in the region 
between approximately 200 end 500 mb has been derived from Voyager data using 
two approaches. The inversion method applied to IRIS data from clear belt 
regions yields 


q * 0.897 t 0.030 

while the combined IRIS-radlo occultat'on method gives 

q s 0.880 t 0.036. 

Fig. 10 and Table 1 compare the helium abundance derived in the present work 
with those determined from previous infrared observations of Jupiter's 
atmosphere by aircraft (Houck et al . 1975; Gautier, et al., 1977) and by 
Pioneer 11 (Orton and Ingersoll, 1976). Abundances are expressed in terms of 
both q and the helium mass fraction. Neglecting the small effect of heavy 
elements, the helium mass fraction is approximately 
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Because of the present uncertainty in the error estimate for the IRIS-radio 
occultation method, we have chosen to use the value derived from the inversion 
method as a basis for comparison in Fig. 10. Within the quoted errors, the 
Voyager IRIS determination agrees with those of the previous infrared 
investigations. The improvement afforded by the present study is actually 
greater than depicted in Fig. 10, since the errors quoted in the earlier works 
omitted estimates of systematic errors, especially those arising from 
inaccurate absorption coefficients. 

The present atmospheric helium abundance of Jupiter may differ from the 
planet's bulk composition, because of helium differentiation during the 
planets's evolution. Differentiation is possible either because helium and 
hydrogen are immiscible over a range of temperatures and pressures relevant to 
Jupiter's interior, or because the metallic-molecular hydrogen transition near 
3 Mbars is first-order, implying a discontinuity in helium abundance across 
the phase boundary. Experiments performed by Streett (1973) on 
hydrogen-helium fluid mixtures at pressures up to 9.3 kbars and temperatures 
from 26 to lOOK have indicated the existence of separated fluid phases. 
Stevenson and Salpeter (1977a, b) have theoretically studied the phase 
behavior of hydrogen-helium fluids at the higher temperatures and pressures 
characteristic of Jupiter's interior and have considered the possible 
conditions for helium differentiation. They conclude that the present 
atmospheric abundance of helium can range from values which are consistent 
with the primordial mix to values significantly less. The uncertainty in the 
relevant thermodynamics makes it difficult to predict when in the evolution of 
the planet differentiation would begin. 

Considerations based on energetics, however, can place constraints on the 
amount of helium differentiation which has occurred within the planet. That 
Jupiter radiates more energy than it absorbs from the Sun is well established 
(see Hanel et al. (1981) and references therein). Two major sources can 
account for this luminosity excess (Flasar, 1973; Salpeter, 1973): first. 
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cooling of the thermil reservoir produced by gravitational energy 
during an earlier quasl-gaseous rtontractlon phase of Jupiter, •'n:d second, 
gravitational separation of hydrogen from hellian during the current epoch. 
Homogeneous evolutionary models of Jupiter, which do not Include helium 
differentiation, predict that the Internal luminosity of Jupiter decays to Its 
currently observed value In 4,5 billion yefrs, the approximate age of the 
solar system (Graboske et al., 1975; Pollack et fil., 1977; Hubbard, 1977; 
Grossman et al., 1980). Significant helium differentiation would Increase the 
decay time by a factor of 5 (Stevenson and Salpeter, 1977b). Homogeneous 
evolutionary models of Saturn, on the other hand, predict luminosity decay 
times of only 2 billion ysars (Pollack et al., 1977; Grossman et al., 1980). 
These results suggest tnat helium differentiation has not, or has only 
recently, begun on Jupiter, but that significant differentiation has already 
occurred within Saturn. The atmospheric helium abundance of Jupiter would 
then be close to its primordial value; Saturn's would be depleted. Grossman 
et al . ( 1980) have estimTited that the depletion on Jupiter, would be less than 
AY = 0.04; Stevenson (1980) has estimated that the depletion In Saturn's 
atmosphere may amount to aY 0.1. 

Whether Jupiter's atmosphere is depleted in helium can. In principle, be 
Independently determined by a comparison with cosmic and solar abundances. 

The present helium content of the universe is composed of a primordial 
component, produced In the big bang an'' the helium subsequently produced in 
stars. The observed large-scale homogeneity and Isotropy of the universe 
suggest that the primordial component la distributed uniformly. The stellar 
component is distributed nonuniformly because the amount of stellar 
nucleosynthesis has varied from place to place. Since nucleosynthesis results 
in a net Increase In helium, a primordial cosmic helium abundance greater than 
that presently measured In Jupiter's atmosphere would provide Indirect 
evidence for differentiation. Greenateln (1980) has reviewed the various 
techniques used to infer the primordial cosmic helium abundance; systematic 
errors in these methods are difficult to assess. One approach is to determine 
the helium content of very old stars. Carney has applied stellar evolution 
calculations to population II subdwarfs and concludes Y = 0.19 ± .04. 

Carney's result, however, Is sensitive to the assumed abundance ratio 0/Fe 
and to the parallax data used to infer the distances of the subdwarfs. 
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Prot* dly thtt mo't raliable daterminatlon of the primordial halium abundance 
has uome from emission line measurements of predominantly extragalactic 
objects. To detormine the primordial helium abundance, Lequeux et al. (1979) 
observed blue compact galaxies and HII regions within irregular galaxies, and 
used earlier measurements of Orion and of the Large and Small Magellanic 
Clouds; French (1980) observed Markerian and Zwicky compact galaxies which are 
spectroscopically similar to HII regions; Talent (1980) observed HII regions 
within trrebular and late-type spiral galaxies. The measurements have 
furnished helium abundances as a function of heavy element abundance. 
Extrapolating the heavy element abundance to zero provides an estimate of the 
primordial helium abundance. Considering the different objects observed, it 
is rather remarkable that the three determinations agree so well (cf. Table 1 
and Fig. 10); Yp Ir 0.22-0.23 with an uncertainty of less than 0.02. Within 
the quoted errors, these abundances are consistent with the current 
determination for Jupiter's atmosphere (Fig. 10). 

A more direct way of determining whether Jupiter's atmosphere is depleted 
in helium is to examine the helium abundance of the Sun, since both bodies 
presumably condensed from the same gas cloud. Observations of helium emission 
lines in the chromosphere and in prominences, and measurements of the relative 
abundances of cosmic rays both yield Y 0.21 with ' typical error of 
approximately 0.04 (cf. Table 1 and Fig. 10). The question arises as to how 
representative these values are of the solar composition. The traditional 
view, from stellar structure calculations, is that the core of the Sun is 
radiative and is surrounded by a thermally convective envelope. No exchange 
of material takes pl'jce between the solar atmosphere and the central region 
where nucleosynthesis ocours. The observable solar atmospheric composition 
therefore represents the initial helium abundance of the Sun, prior to nuclear 
burning, and is directly comparable to -he Jovian nelium abundance. There are 
grounds, however, for questioning the traditional view of internal mixing (see 
below); moreover, Hirshberg (19731 has stressed that large systematic errors 
can arise in the interpretation of the cosmic ray and helium emission line 
data, resulting in uncertainties far in excess of the quoted errors. 

Stellar evolution models, fitting the mass, luminosity, and age of the 
Sun, offer another means of determining its helium abundance. "Standard" 
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models yield initial mass fractions ^ .22-. 24 (Iben, 1969. Bahoall et al.. 
1973; Ulrich and Rood, 1973; Mazzitelli, 1979). A major difficulty with these 
models, however, is their inability to satisfactorily account for the 
anomalously low measured neutrino flux (Davis et al., 1971). Bahoall et al . 
(1973) have fit the low neutrino fluxes with models assuming p low 
heavy-element abundance, and find * 0.115. Since, from the study of 
Stevenson and Salpeter (1977b), the present Jovian atmospheric helium 
abundance places a lower limit on the initial solar helium abundance, this low 
value is precluded by the current work (Fig. 10). If one breaks with the 
standard model and permits sufficiently rapid mixing throughout the solar 
interior, higher values of Y^ are consistent with the observed neutrino flux 
(Iben, 1969). Recently, Schatzmann et al. (I960) have argued that turbulent 
mixing within the radiative core of the Sun can be Induced by the shear 
associated with differential rotation, and have parameterized this mixing in 
terms of a Prandtl mixing-length diffusion term. They are able to fit the age 
and luminosity of the Sun and the neutrino flux with Y^^ = 0.25. The degree of 
internal mixing within the Sun, however, is yet to be examined in a 
dynamically self-consistent fashion and remains an open issue. 

The combined use of interior models with the measured spectrum of global 
oscillations on the surface of the Sun can also furnish an eutimate of Y. 

Iben ’nd Mahaffy (1973) obtained Y^ = 0.19-0.21 in this manner. More recently 
Isaak ^1980), using new observations which resolved the broad 5 minute 
oscillation peak into several lines, concluded Y^ < 0.17, but this limit 
appears to be inconsistent with the estimates of the primordial cosmic helium 
abundance discussed earlier. We conclude that a reliable estima'.e of the 
initial solar heliun abundance remains elusive. The present determination for 
Jupiter may in fact represent the best estimate of the solar value to date. 

To summarize, the estimated errors in the current determination of the 
helium composition in Jupiter's atmosphere, as well as uncertainties in both 
the solar and primordial cosmic abundance of helium, preclude a firm 
conclusion concerning helium differentiation within Jupiter. The results are 
consistent with a uniform mix of hydrogen and helium within Jupiter's 
interior, but a modest amount of depletion in the atmosphere (aY < 0.05) 
cannot be excluded. Comparison with Saturn's helium abundance may prove to be 
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the most valuable diagnostvo tool. An atmospheric helium abundance on Saturn 
lower than that on Jupiter would tend to support the arguments for 
differentiation based on energetics. A helium abundance comparable to that on 
Jupiter would alternatively open . possibility that significant 
differentiation has occurred on neither planet and therefore furnish an upper 
limit on the primordial cosmic helium abunoance. The recent Pioneer 11 
determination of Saturn's atmospheric helium abundance was not conclusive. 

The improved determination anticipated from Voyager IRIS is therefore 
essential for resolving this problem. 
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Table 1. Summary of Helium Abundances 
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Aircraft 
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16-.20 

Houck, et al, 1975 
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.12-. 06 

.2U. 

09-.21 

Gautier, et al. 1977 
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Pioneer 1 1 

.88 

i .06 

.2U. 

09-. 10 

Orton & Ingersoll, 1976 

(c) 

Voyager IRIS 
( Inversion) 

.897 

t .030 

.19 t 

.05 

This investigation 

(d) 

Voyager IRIS 

Radio Occultation 

.880 

t .036 

.21 i 

.06 

This Investigation 


(Egress) 







Solar 
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Helium Emission Lines 

.88 
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.21 ± 

.03 

Iiirayama, 1971 

(f) 

Cosmic Rays 



.20 t 

.04 

Lambert, et al. 1967 
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Standard Interior 

.86 

- .88 

.22 - 

.24 

Iben, 1969: Bancall, 


Models (initial 





et al., 1973; Ulrich & 


abundance) 





Rood, 1973; Mazzitelli, 
1979; Yang, et al. 1979 

(h) 

Low heavy element fit 

.939 


.115 


Bahcall, et al. 1973 

(i) 

Well mixed interior 

.86 


.25 


Schatzmann, et al. 1980 

(J) 

Solar oscillations 

.88 

- .90 

.19 - 

.21 

Iben and Mahaffy, 1976 

(k) 

If II 

>.91 


<.17 


Isaak, 1980 


Primordial 






(1) 

Pop. II subdwarf 

.90 

t .02 

.19 ± 

.04 

Carney, 1979 

(m) 

Extragal. HII emission 

.871 

± .009 

.228 ± 

.014 

Lequeux, et al. 1979 

(n) 

II II II 

.879 

t .009 

.216 ± 

.015 

French, 1980 

(0) 

If II II 

.879 

t .008 

.216 ± 

.013 

Talbot, 1980 
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Figure Captlona 


Figure 1. Pressure-Induced absorption coefficients for molecular hydrogen 
with hydrogen and with helium as the broadening gas. Hie upper part Is 
for T ■ 100 K and the lower for 180 K, 

Figure 2. An Illustration of the sensitivity of the optical depth In the 
Jovian atmosphere to che helium mole fraction. 

Figure 3. Root mean square of the residuals between calculated and mesured 
radiances as a function of the hydrogen mole fractions for one of the 62 
spectra analyzed. 

Figure 4. Histogram of the values of q retrieved from 62 TRIS spectra. 

Figure 5. Estimate of the error in q resulting from a relative error aB/B in 
the H^-He component of the abaorption coefficient. 

Figure 6. ”Warra''(W) and ”cold''(C) thermal profiles derived from the radio 
occultation entry of Voyager 1, compared to a profile retrieved from an 
average of 13 infrared spectra recorded at approximately the same 
atmospheric location (12°S, 63°W). The profiles are derived for a 
hydrogen mole fraction of 0.915. W and C are derived from boundary 
conditions of 155 K at 11.5 mb and of 145 K at 10.8 mb respectively. 

Figure 7. Warm occultation profiles similar that of Figure 4, but derived for 
H 2 mole fraction of 0.86 and 0.9*1. 

Figure 8. Brightness temperature spectra calculated from the ingress 

occultation profile for q = 0.835, q = 0.875 and q s 0.915, compared to 
an average of 13 IRIS spectra measured near the location of the radio 
occultation entry. 


Flgurt 9> Variations In brightnasa temparaturea of tha oaloulatad apaotrum 
raaultlng from variation of tha total absorption coafficiant by ♦SI and 

-51. 


Figura 10. Comparison of halium abundanoa datarmlnations for Jupitar, the 
Sun, and aalrctad aatrophysioal objects. Tha sources for tha various 
values are: (a) Aircraft measurements (Houck at al., 1975: Gautier 

at al., 1977); (b) Pioneer 11 infrared radiometer (Orton and Ingersoll, 
1976); (c) IRIS inversion method (this work); (d) IRIS-radio occultatlon 
method for tha egress point (this work); (a) Solar chromospheric emission 
lines (Hirayama, 1971); (f) Cosmic rays (Lambert, 1967); (g) Standard 
evolution models (Iben, 1969; Bahcall et al., 1973: Ulrich and Rood, 

1973: Mazetelll, 1979; Yang et al., 1979); (h) Low heavy-element model 
fit to neutrino flux (Bahcall et al., 1973); (i) Models with mixing 
throughout the interior (Schatzmann et al., 1980); (J) Iben and Mahaffy, 
1976); (k) Isaak, 1980; (1) Carney, 1979; (m) Lequeux et al,, 1979; (n) 
French, 1980; (o) Talent, 1980. 
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